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The role of large-scale magnetic fields in galaxy NGC 891: can 
magnetic fields help to reduce the local mass-to-light ratio in the 
galactic outskirts? 
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ABSTRACT 

We address the problem of the influence of large-scale magnetic fields on galactic rotation 
for the example of the spiral galaxy NGC 891. Based on its rotation curve and the surface 
density of Hi we determine, in the framework of the global disc model, the surface density 
of matter. Then, based on the surface brightness, we determine the corresponding profile of 
the local mass-to-light ratio. We also model the vertical gradient of azimuthal velocity in 
the quasi-circular-orbit approximation, and compare it with measurements. We discuss what 
factors may influence the rotation of matter in NGC 891 and how this can translate to changes 
in the profile of the local mass-to-light ratio. In particular, we discuss the possible effect 
of magnetic fields on the motion of ionized gas, and, consequently, on the determination of 
the profile of the local mass-to-light ratio. Finally, we put forward the hypothesis that the 
asymmetry in magnetic fields observed in NGC 891 might be responsible for the observed 
anomalous behaviour of the vertical gradient. 

The definitive version is available at 

http : //onlinelibrary. wiley . com/ doi/10 . 1111/ j . 1365-2966 . 2012 . 20447 . x/abstract 
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1 INTRODUCTION 

NGC 891 is a nearby edge-on spiral galaxy. Its very high in- 
clination angle allows the vertical gradient of the azimuthal 
component of the velocity field to be determ ined. Measure- 
ment s include gradients of —15 km s -1 kpc" 1 i Oosterloo et all 
120071) . -18.8 kms" 1 kpc" 1 dKamphuis et al.ll2007h and -17.5 ± 
S.gkms^kpc" 1 for z £ (1.2,4.8} kpc, r 6 (4.02,7.03) kpc 
in the N E quadrant faeald et al.ll2006h . 

In jjatocha et alj|2010h we modelled the vertical gradient for 
the galaxy NGC 891 for the global disc model and quasi-circular- 
orbit approximation, assuming that the motion of gas above the disc 
plane is governed mainly by the gravitation of the column mass 
density distributed in the mid-plane. Our gradient estimate was in 
agreement with the above measurements. This coincidence led us 
to suggest that NGC 891 was a flat, disc-like object, devoid of a 
massive spheroidal component, at least in the region covered by the 
gradient measurements (outside the central bulge). This qualitative 
result concerning the mass distribution is consistent with a more 
recent finding sugg esting that in NGC 89 1 there is no dark matter 
within r — 15 kpc jFraternali et alj|2~01 lh . 

The first aim of this paper is to see to what extent our previous 



E-mail: Lukasz.Bratek@ifj.edu.pl 



gradient prediction is altered by new measurements. To reconstruct 
the column mass den sity we use the rotation curve determined in 
dFraternali et al.l201 lh and the previous me asurement of neutral hy- 
drogen taken from l Oosterloo et al J 120071) . In section [2] we show 
that with the new data our gradient estimate is not significantly al- 
tered. 

Our approac h to mass mode lling is different from the usual 
one presented in dFraternali et alj|2~01 lh . We do not assume a con- 
stant mass-to-light ratio and we do not derive the mass distribution 
based on brightness measurements. First, we obtain the mass distri- 
bution directly from the rotation measurements and the distribution 
of gas. Next, based on the surface brightness measurements, we 
determine local mass-to-light ratios, which then can be drawn as a 
function of galactocentric radial distance. 

It is a feature of the global disc model that the profile of lo- 
cal mass-to-light ratios is sensitive to the accuracy in determining 
the rotation curve. A small variation in the rotational velocity can 
lead to a significant reduction in the mass-to-li ght ratio in the galac- 
tic ou tskirts. On the other hand, it is known (Battan er et al.l 19921 
2008) that magnetic fields can strongly influe nce the rotational fiel d 
of the gas. More importantly, as concluded in (Battaner et al.l 1992h . 
the flatness of rotation curves for the outer parts of discs of spiral 
galaxies could be the result of interstellar magnetic fields. As we 
will see, the inclusion of magnetic fields can help to reduce the 
usually high value of the local mass-to-light ratio at the disc rim. 
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Regarding galaxy NGC 891, there is some peculiarity in the 
behaviour of the vertical gradient of rotation. The gradient vanishes 
in the SW arm. We test our hypothesis that this vanishing could be 
caused by the presence of a relatively intense magnetic field of a 
dozen or so microgauss. 

The final conclusion o f this paper is tha t our r esults are consis- 
tent with the predictions of iFraternali et alj ( 12011!) that there is no 
dark matter within 15 kpc. Furthermore, non-baryonic dark matter 
could be absent also at larger radii, as high local mass-to-light ratios 
(often associated with the presence of dark matter) can be reduced 
to reasonably small values merely by taking into account the inter- 
action of the gas with magnetic fields. We come to this conclusion 
by applying a simplified model. It seems that the non-trivial task of 
rigorously taking into account galactic magnetic fields might turn 
out to be important for our understanding of galactic dynamics. 



2 MODELING OF THE VERTICAL GRADIENT OF 
ROTATION 

We reconstruct the substitute surfac e mass density in the global disc 
model based on the rotation curve ( IFraternali et al.|[201 lh shown in 
Fig[T] and the measur ed surface mass density of neutral hydrogen 
(Ooste rroo et alj |2007) multiplied by a factor of 4/3 to account for 
the presence of helium. The relationship between the rotation curve 
« c (r) and the surface mass density a(r) in disk model, where r is 
the galactoce ntric distance in the disk plane, is given by the follow- 
ing integraf] < Br atek et ai]|2008h 
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Here, x ls an integration variable expressed in the same units as 
r. It is important to remember that when the rotation data are cut 
off, so is the integration region, the result of integration contains 
a boundary term that can differ from that in another integral with 
the de rivative of v\ and given, for example, in dBinnev & Tremaind 
1 1987b . The two integrals are equivalent, provided they are not cut 
off. After the integration has been formally cut off, they give dif- 
ferent results because of the singularities in the integral kernels. 
However, the differen ce is not importan t for radii out to 60 per cent 
of the cutoff radius dBratek et alj|2008f) . This shows that the disc 
model must be used with due care, as the result depends on what 
is assumed about the rotation curve beyond the cutoff radius. The 
integration should not be merely cut off. Moreover, when the in- 
tegration is cut off then the calculated mass density can be lower 
than observed for hydrogen in the vicinity of the cutoff radius. To 
minimize the uncertainty, we find the surfac e density iteratively , 
using a method similar to that described in djalocha et alj|2008f) . 
The method finds a mass density that agrees with the distribution 
of gas for radii greater than the cutoff radius, and the rotation cal- 
culated for this extended mass distribution agrees with the rotation 
observed for radii smaller than the cutoff radius. The results are 



shown in Fig|2] together with the surface density of neutral hydro- 
gen multiplied by a factor of 4/3, separately for the NE and SW 
quadrants. 
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Figure 1. Rotation curve of spiral galaxy N GC 891. The sol id circles denote 
the measurement points taken from i Fraternali et al. 2011), and the line is 
the model rotation curve obtained in the global disc model. 



In the quasi-circular-orbit approximation the rotational veloc- 
ity close to the mid-plane is given by by 
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This relation allows us to determine the vertical gradient of the az- 
imuthal veloci ty, d z v v . In practice it can be calculated in two ways, 
similarly to in (Jaloch a etalj2010h . In method I we calculate v v for 



an array of points (r 



. Then for each Zj we find an average over 



all Ti in the specified range of r. The results are shown in Fig[3] The 
error bars represent a standard deviation. With the use of a linear fit 
the gradient estimate is found to be —19.6 ± 2.9 kms" 1 kpc -1 . In 
method II, gradients are found separately for every r, with the use 
of a linear fit. Next, the mean value of the gradients and the uncer- 
tainty defined as one standard deviation are calculated. In that case 
the gradient estimate is —19.16 ± 1.75 kms -1 kpc -1 . For com- 
parison, we recall the previous result with another rotation curve: 
-lg.g d^kms" 1 kpc" 1 a nd -19. 5± 1.7 kms -1 kpc" 1 , respec- 
tively djalocha et alTl2010h . The result is only slightly changed, 
which is to be expected as the two rotation curves are similar in 
the region of interest. This change is insignificant within uncertain- 
ties. 

It is worthwhile to stress that it is a feature of disc models 
that they give high gradient values. Moreover, our result shows that 
the azimuthal velocity is a linear function of z, and the relatively 
small standard deviation in method II shows that the gradient is 
weakly dependent on the radial variable. These characteristics are 
consistent with measurements. 



1 E and K are complete elliptic integrals of the first and second kind: 

ir/2 tt/2 



3 THE PROFILE OF LOCAL MASS-TO-LIGHT RATIO 



K(ifc) = f — E(fc) = f \/l - k 2 sin 2 <f> dfi ^he P ronle °f me loca l mass-to-light ratio for NGC 891 is shown 

sin 2 cj> J in Fig[4] separately for the NE and SW arms (together with the av- 

eraged profile). To obtain it we used the 3.6/im surface brightness 
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Figure 2. Surface mass density in the global disc model calculated based 
on the rotation curve taken from iFraternali et alj|201 1[) [solid line]; the 
surface brightness in the 3.6^tm filter taken from iFraternali et alj|201 lT) 
shown separately for the NE quadrant [full circles] and for the SW quad- 
rant [ open circles]; the surface density of Hi +He based on ( Ooster loo et alj 
|2007|) [joined full squares]. 
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Figure 4. The profile of the local mass-to-light ratio (without Hi +He). NE 
quadrant [dashed line]; SW quadrant [dotted line]; the mean value [solid 
line]. 
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Figure 3. The azimuthal component of rotation as a function of the vertical 
distance from the mid-plane. The points represent the velocity values aver- 
aged over the interval r € (4.02, 7.03) kpc, and the bars show a standard 
deviation in that interval. 



profile of NGC 891 published in dFraternali etai]|2()I ]]), corrected 
for the inclination of 88° (LEDA), and the surface mass density 
found in the previous section after subtracting the contribution from 
neutral hydrogen and helium (non-luminous in the 3.6/im filter). 
The local mass-to-light ratio is 5 in the central region (outside 1 kpc 
where the rotation and hydrogen data are available), then it falls off 
rapidly to 1 — 2 for r 6 (3,12) kpc and finally rises to a value above 
4 at r — 13.4 kpc. For large radii the surface brightness is known 
only for the NE quadrant. The brightness measurement is not very 
accurate and shows a large spread, making the local mass-to-light 
ratio oscillate rapidly between 1 and 8. 

What conclusions can be drawn? First, the high inclination an- 
gle almost certainly causes the surface brightness profile to be un- 
derestimated (extinction processes can be decisive). Furthermore, 
molecular hydrogen and other gases are present and not accounted 
for - if they were, the local mass-to-light ratio would decrease even 
more. Taking this into account, and also that the ratio is low in most 
parts of the galaxy, we see that there is no need to introduce non- 
baryonic dark matter. The rapid increase in the local light-to-mass 



ratio profile in the outermost region of the rotation curve is also 
not worrying - for large radii in the NE quadrant the average lo- 
cal mass-to-light ratio is above 5, which is not greater than in the 
central region. 

We stress that for large radii the local mass-to-light ratio is 
sensitive even to very small variations in the rotation curve. To 
show this, we change the rotation curve slightly by decreasing it 
in the outer regions, but not by much, so that, except for a single 
point, the profile agrees with measurements within error limits (see 
FigO. Even such a small change can significantly decrease the pro- 
file of local the mass-to-light ratio (see FigHj}. 



250- 



E 150- 




I 1 1 

4 6 



~ 1 1 1 ' 1 ' 1 — 

10 12 14 16 



R[kpc] 



Figure 5. A small modification of rotati on velocities [solid line]. Points 
with error bars were taken from i Fraternali et al. 201 1). 



3.1 The role of magnetic fields 

As we will see below, despite the presence of ionized gas above 
the galactic plane, the influence of magnetic fields with intensities 
of a few microgauss can be ignored in the central region where 
the motion of gas is dominated by gravity. In particular, such small 
magnetic fields should not play any role in the vertical gradient of 
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Figure 6. A change in the profile of the local mass-to-light ratio caused by 
a small modification of the rotation [solid line], compared with the origi- 
nal profile of the local mass-to-light ratio [ dotted line ]. Here, only the NE 
quadrant is considered. 



rotation for NGC 891 in the range r € (4.02, 7.03) kpc. There- 
fore, we expect our gradient model to work well in this region, ac- 
counting satisfactorily for gradient measurements. At larger radii, 
however, the situation is different and including the magnetic fields 
can be essential. 

We have seen that the profile of the local mass-to-light ra- 
tio in the galactic outskirts is strongly dependent on variations in 
the rotation curve. However, the rotation curve in this region is 
determined based on measurements of gaseous mass components 
(mainly hydrogen). If the gas is at least part ially ionized, its mo- 
tion will be influenced b y magnetic fields l lBattaner et al.lfl992l : 
iBattaner & Floridd 1 1995b This effect is especially prominent at 
large radii I Battaner et alj|2008h . For a galaxy with a typical mass 
of about 10 11 Mq, the influence of magnetic fields even at a dis- 
tance of 10 kpc can be significa nt compared with gravitational 
forces dtCutschera & Ja locha 2004). In what follows we will esti- 
mate whether taking this magnetic effect into account can affect 
the local mass-to-light ratio determination. 

Consider an example rotation curve that agrees with the rota- 
tion curve of NGC 89 1 everywhere except at larger radii, where we 
assume the rotation to be slower than observed (see Fig|7J. The re- 
sulting surface mass density is such that the predicted profile of the 
local mass-to-light ratio is a non-increasing function of radius (see 
FigHJ, in contrast to the situation shown in Fig|4]with the original 
rotation curve. Now, let us calculate how strong the magnetic field 
should be to cause such a change in the rotation curve. The starting 
point is the stationary pressure-less and viscous-less NavierStokes 
equation: 



(v o V)v = -V<& 



(V x B) x B. 

Aitp 



(3) 



Here, p is the density of gas. It is assumed that the large-scale mag- 
netic field B is given. Close to the disc, assuming the velocity field 
to be of the form v v (r)e v , and the magnetic field to be of the form 
B v {r)e. v + Bz{r)e. z , the radial component of these equations gives 



YM — 1 

r dr 8np 
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and the other components are identically satisfied. Introducing the 



Figure 7. The measured rotation curve iFraternali et alj201 ll) [circles] and 
rotation velocities lowered so that the local mass-to-light ratio will be a 
non-increasing function of galactocentric distance [ solid line] 
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Figure 8. A non-increasing profile of the local mass-to-light ratio corre- 
sponding to a rotation curve suitably changed for large radii as shown in 
FigCD 



from the gravitational field, we have 

{5v v {r)f _ 1 
r 8wp 



dr 



Here {8v v ) is the contribution to the rotational velocity from the 
magnetic field. First, consider the case of an azimuthal field (B z = 
0). Then the required magnetic field to account for (5v v ) 2 would 
be 



notation V,„ 



+ {8v v ) and denoting by vz, the contribution 



We assumed an exponential falloff of the bulk density p along the 
direction normal to the disc plane: p = poe~ z ^ h , with h = 1 kpc. 
Then po « a/2h, where a is the s urface density of hydrogen HI 
measured in dOosterloo et alj|2007t) (see Fig|9]l. The result of the 
numerical integration of equation l(4j is shown in FigllOl 

It can be seen that a magnetic field of the order of a few mi- 
crogauss suffices to produce the required difference in velocities. A 
similar calculation can be repeated assuming a purely vertical mag- 
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Figure 9. Bulk d ensity of Hi in t he disk plane obtained based on the surface 
density given in I Ooster loo et alj|2007l) . 




Figure 10. Magnetic field needed to lower the rotation so that the profile of 
the local mass-to-light ratio is a non-increasing function of the galactocen- 
tric distance: the vertical component B z of magnetic field [dotted line], and 
the azimuthal component B,p of the magnetic field [solid line]. 



netic field with B z being a function of the radius only. In that case, 
equation ([3} simplifies to 



(X \ 2 r 9 n2 



with the solution 



(B z (ri)) 2 + 87r 



p(0(*v(e)) a £de- 



(5) 



(6) 



In equation l[6ll we assume the same density p as in equation {4j. 
The result of the numerical integration is presented in Fig[l0] Sim- 
ilarly to in the previous case, the field is of the order of several 
microgauss. In the galaxy NGC 891 there is indeed a large-scale 
magnetic field presen t of the order of a few microgauss (both in 
the disc and above it) ( Alle n & SukumaJ 19911 : iDumke et al J 19951 : 
iHummel etaljl99lh . 

With the magnetic field as shown in Fig[l0] the mass-to-light 
ratio profile remains constant, but weaker fields can reduce this ra- 
tio. The plots in Fig[l0]are shown for illustrative purposes only. 
They allow estimation of the strength of the magnetic field required 
to reduce the mass-to-light ratio. 



The situation considered above serves as a toy model, illustrat- 
ing the fact that the presence of magnetic fields with an intensity of 
a few microgauss cannot be simply ignored, as they have a sig- 
nificant impact on the motion of (even only partially) ionized gas. 
Not accounting for this effect may lead to an overestimation of the 
profile of the local mass-to-light ratio. 



4 THE VERTICAL GRADIENT OF ROTATION AND 
MAGNETIC FIELD 

We recently applied our gradient estimation method to six galax- 
ies with available measure ments of the vertical gradient of rotation 
dJalocha et al.ll2010l . l201 11) . We found our predictions of the gradi- 
ent value to agree within error bars with measurements. This agree- 
ment was reached without taking into account the influence of mag- 
netic fields. This is so because the gradient was measured only at 
small distances from the galactic centres (eg. 4— 7kpc for NGC 891 
or 3— 8kpc for Milky Way). As will become evident later, and as 
follow s from a simple estimation made in (Kutschera & Jalochal 
2004), magnetic fields with a typical intensity of several micro- 
gauss are too weak, compared with the gravitational influence, to 
have a noticeable effect. For gradient measurements at low altitudes 
above the disc plane (as for the Milky Way), the influence of mag- 
netic fields can be neglected with greater certainty when the gas 
density decreases with the altitude above the galactic plane. In Sec- 
tion [2] we had such an agreement with gradient measurements for 
galaxy NGC 891 in the NE quadrant. 

There can be stronger magnetic fi elds present in galax ies, of 
the order of a dozen or so microgauss iBattaner et ai]|2008h . Such 
strong fields can be decisive for the motion of ionized gas, in par- 
ticular of gas of low density. A gradient is not observed in the SW 
arm of NGC 891. Could a magnetic field be responsible for this 
anomaly, reducing the gradient to zero? To answer this question 
we again employ the NavierStokes equation l[3). For simplicity, we 
can assume a purely azimuthal magnetic field that depends on the 
variable z only. Then: 



vl(r,z) = -r^-$(r,z) 



+ 



where "&{r, z) is the rotation velocity observed above the disc, and 
v^(r, z) is the rotation velocity one would expect in the absence of 
a magnetic field. A constant gradient value is observed above the 
mid-plane, starting from z. a = 1 kpc. Hence, S|(r, z) = v^(r, z a ) 
for z > z a . This gives us 

B\ = B\(z a ) + Anp (vl(r, z a ) - v%{r, z)) . (7) 

The density p in this equation is shown in Fig[TT] We assumed an 
exponential falloff of the density along the direction normal to the 
disc plane, with h = 1 kpc. In fact, the obtained is also a func- 
tion of r; however, it is weakly dependent on r compared with the 
changes associated with the variable z. Hence, the result for B^ 
can still be regarded as consistent with the earlier assumption that 
d r B,i, — 0. The behaviour of B^ is shown in Fig|12l 

It can be seen that to reduce the gradient of rotation in 
the region r 6 (4, 7) kpc, a magnetic field with the intensity 
of a dozen or so microgauss would suffice. This means that we 
need fields larger than t ypical for spiral galaxies. T he results of 
lAUen & Sukumarl Jl99lh and IHummel et all dl99ll) suggest that 
particularly strong and ordered fields might indeed be locally 
present in the SW arm of galaxy NGC 891. Similarly to in the pre- 
vious section, we stress that the field values are shown for the pur- 
pose of illustration only. They show the order of magnitude of the 
magnetic field required to cancel the gradient. With this in mind, 
we can conclude that a field of a dozen or so microgauss can affect 
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Figure 11. The bulk density of HI as a function of vertical distance from 
the disk plane. 
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Figure 12. Magnetic field required to cancel the vertical gradient shown at 
r = 4.0 kpc [solid line] and at r = 7.0 kpc [dotted line]. Calculations of 
the field assume the gas density as in Fig ll II 



the gradient and, if there were such a field in the SW arm of NGC 
891, it could be used as a possible explanation for the anomalous 
gradient. There is indeed a certain field asymmetry observed: 'in 
the south-western part of NGC 891 the polarization percentage is 
between 20 and 30 p er cent, implying very well ordered B -field' 
faummeletai1ll99lh . If the asymmetry applies also to the field 
intensity, then we would have a strong premiss for attributing the 
gradient asymmetry to this magnetic effect. 



5 SUMMARY 

From the above it follows that galaxy NGC 891 is a disc-like ob- 
ject without a dominant massive non-baryonic dark matter halo. 
Two facts support this statement. First, the vertical gradient values 
of the azimuthal velocity calculated in the framework of the disc 
model are high, consistent with the measurements. Second, the ob- 
tained profile of the local mass-to-light ratio is low. Although the 
ratio rises for larger radii, its behaviour, as we have shown, is sen- 
sitive to small variations in rotational velocities. To illustrate this, 
we reduced the rotation by a small amount (allowed by measure- 
ment errors). Moreover, the ratio profile may be overestimated by 



the omission of factors such as the presence of gases other than HI 
or the influence of magnetic fields. 

From the above it follows that galaxy NGC 891 is a disc-like 
object without a dominant massive non-baryonic dark matter halo. 
Two facts support this statement. First, the vertical gradient values 
of the azimuthal velocity calculated in the framework of the disc 
model are high, consistent with the measurements. Second, the ob- 
tained profile of the local mass-to-light ratio is low. Although the 
ratio rises for larger radii, its behaviour, as we have shown, is sen- 
sitive to small variations in rotational velocities. To illustrate this, 
we reduced the rotation by a small amount (allowed by measure- 
ment errors). Moreover, the ratio profile may be overestimated by 
the omission of factors such as the presence of gases other than HI 
or the influence of magnetic fields. 

Our gradient modelling accounts well for measurements in 
the NE quadrant. Regarding the SW quadrant, the absence of gra- 
dient could be associated with the presence of stronger magnetic 
fields, with intensities of a dozen or so microgauss. If such asym- 
metry occurs for a galaxy, one can hypothesize that it is caused 
by some asymmetry in the intensity or confi guration of the large- 
scale magnetic fiel d. Based on the research in dHummel et al J 199 it 
I Allen & Sukumarlll99ll) it seems there is such an asymmetry in 
NGC 891. We stress that our considerations concerning magnetic 
fields are illustrative in character. The influence of magnetic fields 
on the motion of matter depends on the field intensity, its configu- 
ration, and on the degree of gas ionization and the gas density. On 
the one hand, accurate determination of these factors is not possi- 
ble; on the other hand, the influence of magnetic fields on galactic 
rotation cannot be ignored. 

In galaxy NGC 891 not only are there are no visible effects of 
dark matter out to 15 kpc ^Fraternali et al.l201 lh . but most probably 
there is also no massive non-baryonic dark matter halo, although 
we cannot rule out its presence completely. 
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